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I. INTRODUCTION 


This thesis investigates the production of nitric oxide 
(NO) by solar soft X-rays in the E-region of the ionosphere 
over the solar cycle. The variation between extreme 
ultraviolet (EUV) flux and X-ray flux as a function of solar 
cycle is investigated as the principal factor influencing NO 
production in the E-region. This research was prompted by 
the hypothesis of Barth et al. (1988) that nitric oxide 
variations are produced by variations in the solar soft 
X-ray flux between 20 and 100 A. The main topic will be an 
analysis of the production of nitric oxide at 110 km as 
measured by the Solar Mesosphere Explorer (SME) satellite. 
These measurements are compared to the results of a 
one-dimensional diffusive photochemical model developed by 
Cleary (1985) at the University of Colorado and modified by 
Siskind (1988) to allow for large variations in the solar 


X-ray flux between 18.62 and 50 A. 


A. THESIS GOALS 

This thesis is a theoretical and observational study of 
the production of nitric oxide in the E-region of the 
ionosphere based on variations between EUV and X-ray flux 
over solar cycle. An empirical expression relating X-ray 
and EUV flux is sought in order to better model _ the 


production and loss of NO in the lower thermosphere. 


B. THESIS OUTLINE 

The thesis is divided into seven chapters. Chapter II 
presents background information about the composition and 
structure of the ionosphere, and describes some techniques 
for measuring the density profile of various constituents. 
Atmospheric chemistry is reviewed in Chapter III, with the 
emphasis on photochemistry. The continuity equation and 
physical transport mechanisms are discussed, as well as the 
method used by the model in solving the continuity equation. 
The chemistry of nitric oxide is focused on in section C, 
and the concept of X-ray scaling is introduced. Chapter IV 
discusses photochemical modeling and describes the basic 
Structure and operation of the photochemical model used in 
this thesis. Chapter V proposes an empirical relationship 
between the X-ray scaling factor (FRACXR) and 10.7 cm flux 
over solar cycle, with results of comparisons between model 
data and SME satellite measurements presented in Chapter VI. 
Chapter VII summarizes the thesis and presents suggestions 


for further research. 


II. BACKGROUND 


Nitric oxide makes up only 0.5 parts per billion of the 
Earth’s total atmosphere (NOAA, 1976), yet its role as a 
catalyst in the destruction of stratospheric ozone makes it 
one of the most important trace constituents. The 
1onization of NO by geocoronal Lyman-a is largely 
responsible for maintaining the nighttime E-region 
ionization at about 100 km. Additionally, because nitric 
oxide is very reactive, it plays an important role in the 
daytime ionospheric character. Because this thesis concen- 
trates on ionospheric NO production, a brief description of 


the ionosphere is presented. 


A. THE IONOSPHERE 

Above about 60 km, there are sufficient ions’ and 
electrons to create a region which affects the propagation 
of radiowaves. TitDemregqronmrs Called the 19onosphere. The 
ionosphere is neutral overall. It is composed mostly of 
neutral molecular and atomic gasses which have’ concen- 
trations that are several orders of magnitude greater than 
the ions and electrons. The variation of electron density 


with altitude leads to three natural subdivisions within the 


ionosphere. These three regions are called the D, E, and F 
regions. The F-region is further divided into’ two 
Subregions called the Fl- and F2-regions. This character- 


istic "layering" of the ionosphere is due to the variation 
of atmospheric composition with altitude and the wide range 
of ionizing wavelengths associated with each neutral 
Species. Figure 2.1 shows the different regions of the 


ionosphere. 
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Figure 2.1. Regions of the Ionosphere. 


1. The D-Region 
The D-region is the lowest lying ionospheric region 
and is therefore produced by ionizing radiation at those 
frequencies which penetrate deepest into the atmosphere. It 
extends from about 50 to 90 km in altitude. The major 
source of ionization in the D-region is photoionization of 
NO by the Lyman-« 1215.7 A line. 
2. The E-Region 
The E-region lies between about 90 and 130 km, and 
has a daytime electron density about two orders of magnitude 
greater than that of the D-region. Ionization in this 
region iS primarily due to photoionization of molecular 
oxygen and nitrogen by EUV and soft X-rays. 
3. The F-Region 
all 


In the F-region, above 130 km, O, O and N 


i hi 2 
photoionize, with the primary source above 150 km being 
atomic oxygen. Within the region, the electron density is 
fairly uniform until about 160 km, where it then rises to a 
daytime maximum of approximately 10° /em? at 220-450 km. 
Because of this characteristic profile, the F-region is 
further divided into the Fl- and F2-regions lying at 130-160 
km and 160-600 km respectively. In the Fl-region, ion 
chemistry dominates the characterization of the ionosphere 


whereas in the F2-region physical transport (1i.e., dif- 


fusion) is the most important process. 


B. THE THERMOSPHERE 

The Earth’s neutral atmosphere is divided into regions 
primarily on the basis of temperature gradirents. Typicaumia 
at middle latitudes, the temperature falls with increasing 
altitude until reaching a local minimum of around 200 K at 
about 10 km. Here the atmosphere undergoes a temperature 
inversion where the temperature rises with increasing 
altitude until a local maximum of 290 K is reached at about 
50 km. The atmospheric temperature then again falls with 
increasing altitude until an absolute minimum of about 160 K 
is reached at 90 km. This minimum is termed the mesopause 
and marks the beginning of the thermosphere. In the themmge 
sphere, temperature again increases with altitude, 
approaching a maximum value in the region of 1000-2000 K at 
a height of several hundred kilometers. Above the 


thermosphere, the atmosphere is essential lyvercotnemmnaas 


C. REMOTE SENSING TECHNIQUES 

Because the ionosphere is a plasma, it will affect the 
propagation of radio waves. It is for this reason that 
radio sounding methods were the First experimental 
technigues used for studying the ionosphere (see e.g., 
Bauer, 1973). Other techniques can be used to obtain 
experimental information about the chemistry and chemical 
processes of the atmosphere. Among the most important of 
these are ground-based observations of airglow emissions and 


the use of satellite-borne photometers and spectrometers. 


1. Remote Sensing of Chemical Processes 
a. Ground-based Airglow Observations 
The term airglow is used to describe non-auroral 

light emissions from the upper atmosphere. These emissions 
are extremely weak and therefore require very sensitive 
equipment to be detected. A common instrument used for 
airglow measurements is the scanning monochromator which 
detects light in a very narrow bandwidth. A major 
difficulty encountered in airglow detection is the isolation 
of the emission from the strong background of scattered 
sunlight. (see e.g., McEwan and Phillips, 1975) 

b. Satellite Observations 

Satellite airglow observations offer an 
advantage over ground observations in that there is almost 
no scattered sunlight to obscure the actual airglow 
emissions. An additional advantage lies in the geometry of 
Satellite limb viewing (see e.g., Liou, 1980). This makes 
possible the calculation of the column density of an 
emitting species. Through an "onion-skin" inversion method 
the number density of the emitting species can then be 
calculated (see e.g., Cleary, 1985; Hickman et al., 1977). 
2. Radio Sounding of the Ionosphere 

Electron density profiles of the ionosphere can be 

measured uSing a variable frequency radar which measures the 


time delay between the transmission of an RF pulse and its 


echo from a reflecting level within the ionosphere (see 
e.g., Bauer,1973). Such a device is called an ionosonde. 
The requirement for reflection of an electromagnetic 
wave from an ionized medium is that the magnitude of the 
complex propagation vector, kK, be imaginary. This require- 


ment is met when 
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The value of w for which equation (2.1) is identically one 
is Wy the plasma frequency. The index of refraction, n, in 


the ionized medium is given by 
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The requirement for reflection of an electromagnetic wave at 
normal incidence is n= 0. Thus the maximum frequency for 


normal reflection is given by 


f° = 80.6{e] (2.3) 


where £ is measured in Hz and [e] is the electron density in 
electrons/m>. Thus by varying the frequency it is possible 
to map out the electron density profile. However, 
examination of Figure 2.1 shows that the electron density 


above about 250 km cannot be measured by ground-based 


ionosondes because the frequencies needed to probe this 
region will already have been reflected at a lower altitude. 

In order to obtain electron density profiles for the 
region above 250 km, ionosondes have been placed. on 
Satellites. Top-side sounding is a very valuable technique 
for obtaining electron density profiles in the upper 
F2-region. Unfortunately, most of the measurements taken 
from top-Side ionosondes do not provide full profiles down 
to the peak altitude so that it is difficult to merge this 


data with that of ground-based ionosondes. 


IIIT. CHEMICAL REACTIONS IN THE THERMOSPHERE 


A great variety of atoms and molecules exist in the 
thermosphere. As will be discussed later, much of this 
molecular soup is subject to reaction with solar radiation. 
The resulting photolytically generated species undergo many 
different chemical processes with the surrounding 
atmosphere. This chapter will examine these processes and 
introduce the method used for analyzing the rapidly changing 


composition and structure of the thermosphere. 


A. PRINCIPLES OF PHOTOCHEMISTRY 
Photochemistry 1s concerned with chemical reactions that 
are initiated by the absorption of solar radiation. The 


absorption process is written as 
* 
A+ hveraA, (3.9 


where Ae 1s an excited energy state of molecule A. The 
excited state could be a rotational or vibrational level, an 
electronic state, or it could represent the decomposition of 
the molecule into two or more fragments. This latter 
process 1s called photodissociation or photolysis (ites 
products are all neutral atoms or molecules, and is called 
photoionization if the products include an ion and) omemias 


more electrons. 
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Once excited molecules or atoms are produced, they 
undergo secondary chemical reactions, forming species that 
would not exist in a normal state of local thermal 
equilibrium (LTE). Therefore these photochemical processes 
are non-LTE mechanisms and must be quantitatively analyzed 
by studying the rate at which reaction occurs (see e.g., 
Chamberlain,1978). 

1. Photodissociation 

The process by which the absorption of a photon 
leads to disruption of the absorbing molecule is called 
photodissociation or photolysis. The atoms produced in this 
manner are highly reactive. Photolysis processes which take 


place in the photochemical model used for this thesis are: 


NO + hv = 0 + N(°S), ieee 
4 4 
N5 + hv > N( ' S) + N(S), C353) 
and 
2 2 
N, + hv > N('D) + N(D). (3. 4) 


Notice that the last two equations involve the 
photodissociation of molecular nitrogen into two different 
electronic states of elemental nitrogen. The extent to 
which these two reactions occur is given by their branching 
ratios. For the above two reactions, the branching ratio 


for each has been determined to be 0.5, indicating that half 
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of the reactions which occur wi Pi produce n( 45s) and that the 
other half will result in N(“D) production (see e.g., 
Heiklen, 1976). 

a. Photodissociation Rates 


There are three factors influencing the rate of 


photodissocidael om. : molecular concentration, photon fiia7e 
and quantum yield. The most obvious factor must be the 
amount, or concentration, of molecule A. Concentrations OF 


a substance per unit volume are denoted by enclosing the 
substance in square brackets. [A] is the concentration of 
substance A per unit volume. 

Another factor influencing the rate is the 
number of photons available to cause photodissociation. The 
flux of solar radiation outside the atmosphere is simply a 
function of the brightness of the Sun and the distance from 
the Sun. However, as the radiation penetrates the Earth's 
atmosphere the flux of photons is attenuated due _ to 
absorption by atmospheric gasses. At high altitudes the 
atmospheric density is low and therefore little absorption 
takes place. As the altitude decreases the density 
increases, and absorption becomes important. 

The efficiency, or quantum yield of the  ,eameiwag 
given by equation (3.1) is defined as the number of A” 
molecules produced for each absorbed photon. Quantum yield 
varies with wavelength and is zero longward of the 


lonization wavelength. 


ie 


A common way of expressing the photon flux at a 
particular altitude z is in terms of the exoatmospheric flux 


adetne opumcala depth. Optical depth, t, 1s defined as 


A 
il 


Doe, (i) | [ildz’, (3.5) 
: Z 
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where o. is the absorption cross section for the en 
constituent at a specific wavelength, and [i] is the 
concentration of that constituent. The photon flux is then 


given by 


-t-sec(x) 


$(z) = de (3.6) 


phere) @(Z) 1s Che photon flux at altitude z, $, is the exo- 
atmospheric photon flux, t is optical depth, and x is the 
solar zenith angle. Note that t-sec(x) is the optical 
pathlength through the atmosphere. 
2. Photoionization 
When the products of equation (3.1) include molec- 
ular or atomic ions and one or more electrons, the process 


is called photoionization. Schematically this is written 


+ * 
A+hv®P®A +e. (32) 


In this process, the energy carried in by the photon is 
converted into ionization energy plus kinetic energy of the 


electron. TRescemencrgetac electrons, are called photo- 
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electrons, oa and enter into several different energy 
transfer reactions with neutral gasses. 

The major source of ionizing radiation in the 
thermosphere is solar radiation. The only wavelengths which 
penetrate to the E-layer are extreme ultraviolet (EUV), 
X-rays, and Lyman-ca. 

a. Photoionization Rates 

The rate at which photoionization takes place 
for each constituent is described by its ionization rate. 


For a given constituent, the lonization rate is given by 


ay(z) = | oj (A804) ad, (3.8) 


where o is the ionization cross section and $. is the solar 


Flux. at valtatudesz. (Torr et al., 1979) 


B. DIFFUSION AND THE CONTINUITY EQUATION 

It is possible to construct an equation showing the time 
variation of a particular constituent’s concentration once 
the rate coefficients for the production and loss mechanisms 
are known. This equation is called the continuity equation 
and contains terms relating the chemical creation and 
destruction of the constituent, as well as terms describing 
the flux due to physical transport processes. There are 
four types of physical transport mechanisms: molecular 


diffusion, ambipolar diffusion, eddy diffusion, and thexiad 


dtitius 1On. 
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1. Molecular and Ambipolar Diffusion 
When the concentration gradient of a constituent in 
the atmosphere differs from its equilibrium value, the 
resulting motion is called molecular diffusion. For an 


isothermal atmosphere, this can be written 


= afi] [i] 
Ds (2s + =| (3.9) 





ee 


where $5 is the flux, D, is the molecular diffusion coef- 


th constituent. 


ficient, and H; is the scale height for the i 

Diffusion which results from a concentration 
gradient of a charged constituent in a partially ionized 
gas, or plasma, in which the constituent is influenced by 
the electrostatic forces of the plasma as a whole is called 
ambepolar Gitfusion. Although electrons in the thermosphere 
have a greater diffusion coefficient than the ions, 
electrostatic attraction between electrons and positive ions 
prevents any large scale separation. Thus the plasma 
diffuses as a whole, with ambipolar diffusion coefficient 
D.: (see e.g., McEwan and Phillips, 1975) 

2. Eddy Diffusion 

Eddy diffusion is the result of turbulent mixing in 
the presence of a composition gradient. Whereas molecular 
diffusion results from the average velocity of the 


constituent being different from that of others, eddy 


diffusion occurs whether or not there is any difference in 
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average molecular velocity between the constituents. Eddy 
diffusion can be expressed similarly to molecular diffusion 
by an equation analogous to (3.8). The coefficient of 
molecular diffusion D; is replaced by the eddy diffusion 
coefficient K. 

3. Thermal Diffusion 


Thermal diffusion is due to large temperature 


gradients in a gas mixture. It can be properly accounted 
for only in a rather complex manner. However, since the 
thermal diffusion factor a. is proportional to- the 


1 


difference in molecular weights of the constituent gasses, 
thermal diffusion is important in the atmosphere only for 
light gasses such as hydrogen and helium. Because the model 
used in this thesis 1S primarily concerned with the 
chemistry of  mnieGrae weexide, thermal diffusion can be 
neglected. 
4. The Continuity Equation 

Physical transport and chemical processes act on 

each constituent of the atmosphere to produce a variation in 


concentration with time. The equation which governs these 


processes 1s the mass continuity segquacueen 


d[1] 
dt 
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+ [Cil-v, | =P, - bli). (3.10) 
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Again, {1i] ais’ the concentration at time t and 
altitude 2Z, P. and L. [1] are chemical production and loss 
rates, and V; is the mean vertical velocity due to transport 
processes. Since [i]-v, is the vertical flux ei, the second 
term in equation (3.9) may be written 

3 ori 
-(til-v,} = >>. (3) 
For a neutral constituent with low concentration, $. may be 


written as the sum of molecular and eddy diffusion. Thus 
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where the middle terms are due to the dependence of scale 
height on temperature gradient, D, and K are the molecular 
and eddy diffusion coefficients respectively, T is the 
absolute temperature, H; the scale height for the jen 
constituent, and Cae the scale height for a constituent 
having the mean molecular mass. (Colgrove et al., 1966) 
a. Time Constants 

The continuity equation includes chemical and 
dynamical processes. The dominance of one process over 
another can be determined by examining the relative 
Magnitude of the characteristic times, or time constants, of 


the competing processes (Goody and Walker, 1972). If two 


processes are competing and one acts more rapidly than the 


so 


other, the faster one will be more effective. The dominant 
process will be the one with the shorter time constant (see 
e.g., Wallace and Hobbs, 1977). 

By performing dimensional analysis on _ the 
continuity equation, the time constants for chemistry, Tor 
eddy diffusion, Tres and molecular diffusion, Ty, can be 


approximated (Cleary, 1985). The distribution of [ij is 


< tT, and tt. <i 


determined principally by chemistry if To K c Dp! 


and by dynamics if TT. < T, and Ty < Te: 
A comparison of the lifetimes for molecular 


diffusion, eddy diffusion, and chemical destruction was done 


by Cleary (1985) for the different species used in that 


model. The comparison showed that dynamical processes are 
only important for two species, NO and n(4s). For the 
remaining species (N(“D), NO*, 0,*, N,*, 0*, and 0(7pD)*), 


chemical processes are sufficiently fast that transport may 
be neglected (Cleary, 1985). Ambipolar diffusion for o* was 
later added to the model. 
5. Solutions to the Continuity Equation 
For constituents which have sufficiently short 
chemical time constants, the dynamical processes’) of 
diffusion may be neglected. Thus the flux term in equation 


(3.10) vanishes and the continuity equation reduces to 


i (3.10a) 





O , and o(“p)*. 
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For those constituents which have short dynamical 
time constants (i.e., NO and N(?s)), the flux term in the 
continuity equation must be retained. When the expression 
for @ given by equation (3.12) is substituted into equation 
(3.10) and differentiated, a parabolic differential equation 


results. This equation can be written 





Doss ; 
d{iJ _ 9 {1] d{i]} ’ 
at = A; + 2 a + Cc, {il + E.- 3.13) 
dZ 
The coefficients A., B., C., and E. in equation 


(3.13) are given in Table I. Equations (3.10a) and (3.13) 
cannot be solved analytically and must be solved by 
numerical methods. In the model used for this thesis, the 
method is the Crank-Nicolson finite differencing technique. 


(Cleary, 1985 App. A) 


C. CHEMICAL REACTIONS IN THE LOWER THERMOSPHERE 

Many types of chemical reactions are continually taking 
place in the atmosphere. Products of these reactions include 
photolytically generated atoms, radicals, excited molecules 
and ions. As was seen in previous sections, photoionization 
and photodissociation are the primary photochemical 
reactions. Reactions involving the products of these two 
types of reactions with themselves and with the _ stable 
atmospheric gasses determine in great part the 


characteristics of the atmosphere. 


JL, 


TABLE I 


COEFFICIENTS TO THE PARABOLIC DIFFERENTIAL EQUATION 





The appendix summaries the types of chemical processes 
which are important in the execution of the one dimensional 
model used in this thesis for the calculation of nitric 
oxide in the thermosphere. Of course, many other reactions 
also occur, but it is the chemistry of nitric oxide which is 


germane to this investigation. 
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1. Nitric Oxide Chemistry 

Direct sources of nitric oxide in the E-region are 
the reactions of atomic nitrogen (ground state n( 4s) and 
excited states (7D) and N(“P)) with molecular oxygen. 
These three reactions (R28, R29, and R39 of Appendix) do not 
equally contribute to the production of NO. In steady 
state, the NO production rate of each reaction is found by 
multiplying the rate coefficient by the product of the 
concentrations of the reactants. Figure 3.1 shows some 
typical daytime nitrogen density profiles. Figure 3.2 shows 
the corresponding NO production rates from the three 
different atomic nitrogen states, as well as the _ rate 
coefficient (k) for each reaction. It can be seen from 
Figure 3.2 that the principal source of NO at 110 km is 


N(“D) via the reaction 


N(7D) + 0. » NO +0 (3.14) 


Z 


The small contribution from N(“S) is due to the temperature 
dependence of the rate coefficient, and the small 
contribution from N(“P) is because of its low concentration. 

Sources o£ n( 7D) are from the dissociative 
recombination of N,” and NO” (R19 and R27), the dissociative 
excitation by photoelectrons of N. (R12), the photo- 
dissociation of N. (R15), and the ion-atom interchange 
between No” and O (R17). In a method analogous to that used 


for NO production, it can be shown that the two principal 
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Figure 3.1. Atomic Nitrogen Densities. 
sources of N(“D) are via R27 and R12. These equations are 


stated here for ease of reference 


Not + e»> O + N(“D) R27 


+ 


* 2 71 
N, +e 7% N('D) + N('D) + e R12 


The two loss processes for nitric oxide in the E-region are 
via the bimolecular reaction with n(4s), (R33), and the 


charge exchange with Ora C507a) 


n(4s) + NO 9 N, +0 R33 


Ong + NO » NO’ + ore R22 
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Figure 3.2. NO Production Rates 


Thus, while N(*D) is a source of nitric oxide at 100 km, 
N(4s) is both a source and sink. An analysis of these 
production and loss processes as they pertain to the 


photochemical model is presented in Chapter V. 
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IV. THE PHOTOCHEMICAL MODEL 


A. PHOTOCHEMICAL MODELING 

Photochemical modeling is a useful tool for gaining an 
understanding of the minor constituents in the thermosphere. 
By comparing a theoretical model with observations of some 
constituents for which data are available, one can determine 
the validity of the model and better understand physical and 
chemical mechanisms occurring in the thermosphere. 

"Photochemical models are sensitive to various input 
parameters such as solar fluxes, photoabsorption 
cross-sections, chemical reaction schemes and reaction rate 
constants, and transport parameters (e.g., eddy diffusion 
coefficient) as well as to the method of numerical analysis 
and boundary conditions." (Shimazaki, 1985). Additionally, 
models may be time dependent or time independent. Time 
independent models are usually used to arrive at a 
Steady-state condition. A steady-state model does not 
represent the actual state for the constituents whose 
concentrations change as a function of time. Causes for 
these changes would be solar flux, geomagnetic Ap index, and 
any temperature dependent reaction rate since the temp- 
erature is also time dependent. Time dependent models 


attempt to account for these temporal variations. 
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B. THE TIME DEPENDENT PHOTOCHEMICAL MODEL 


1. Numerical Technique 
In Chapter III it was shown that the three phenomena 
which must be understood in order to model an atmospheric 
constituent properly are chemistry, eddy diffusion, and 
molecular diffusion. Additionally, it was shown that the 
continuity equation reduces to equation (3.10a) for species 
which do not require transport processes to be considered. 
For NO and N(7s) transport processes must be included and 
the continuity equation becomes a parabolic differential 
equation (equation (3.13)). 
Equation (3.10a) can be solved by the following 


semi-implicit finite difference equation 


iieldt! eed eae Te 
t 


BE t a) 


t+1 ’ 


which is solved for [iJ,i4) where [1] 
th 


ta] is the density of 


the i constituent at time interval t+l and is a function 
of [il],. 

Equation eS) All Si) is first order in the time 
derivative and second order in the spatial derivative and 
therefore requires three boundary conditions: one in time 
and two in space. Cleary (1985) has shown and Siskind 
(1988) has verified that the spatial boundary conditions are 


satisfied by specifying conditions at the upper and lower 


altitude limits of the model. Several different spatial 


Z5 


conditions can be used to implement these boundary 
conditions. The ones used by the model require that there 
be zero flux of both NO and N( 7s) at the upper boundary, 
and at the lower boundary the condition for NO is that there 
is a local minimum in the density profile, while for n(45s) 
the requirement is that it be in a state of photochemical 
equilibrium. 

When equation (3.13) 1S cast into its finite 
difference form using the Crank-Nicolson method, a series of 
coupled algebraic equations are produced which can _ be 
expressed as a tri-diagonal matrix equation. The density 
profile for each time increment is obtained using the Thomas 
algorithm of matrix inversion. (Von Rosenberg, 1969; Cleary, 
1985) 

2. Inputs to the Model 

The Cleary time dependent photochemical model is run 
in three distinct steps. The first step uses the MSIS-83 
program (Hedin, 1983) to calculate the neutral atmosphere. 
The second step uses a modification of the NRL photoelectron 
code PEGFAC (NRL Memorandum Report 5004, 1982) to calculate 
photoelectron impact rates. Siskind (1988) modified the 
PEGFAC code to extend the energy grid so as to include soft 
X-rays up to 666 ev (18.62 A). The existing energy grid was 
also updated to reflect more recent data. The third step is 


the actual one dimensional model. The model calculates the 
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time development of each species which enters into the 
chemistry of nitric oxide. 
a. Neutral Atmosphere 

The first step of the photochemical model uses 
the MSIS-83 (Hedin, 1983) code to calculate the neutral 
atmosphere for each model "hour" of every day for which the 
model is run. Results are stored as 183 element column 
vectors. Each "hour" five of these vectors are generated. 
The first contains information on the geophysical parameters 
used. These are (1) the exoatmospheric temperature, (2) 
solar zenith angle, (3) the F10.7 value for the previous 
day, (4) the 90 day average F10.7 value, (5) geomagnetic Ap 
index, (6) latitude, (7) longitude, and (8) the X-ray 
scaling factor, FRACXR. The second vector contains the 
temperature profile for altitudes from 86 to 450 km in 2 km 
increments (hence 183 elements). The third fourth and fifth 
column vectors contain the density profiles of 0O, 05, and 
No, respectively. Have cach Nour’ ~1is a 183 x 5 array. 
Each day produces 24 such arrays. 

b. Photoelectron Impact Rates 

The calculation of photoelectron impact rates 
was done using the computer program PEGFAC (NRL Memorandum 
Report 5004, 1982). Recently, Siskind (1988) revised the 
PEGFAC code to account more properly for the higher energy 
photoelectrons produced by soft X-rays. He expanded the 


Original 37 EUV intervals from 50-1050 A (Torr et al., 1979) 


Jy | 


to include an additional 20 intervals which cover the soft 
X-ray flux from 18-50 A. The fluxes used were measured in 
one of Hinteregger’s rockets and are available from the 
National Space Science Center ina file labelled SC#21REFW 
(Siskind, 1988). 

EUV fluxes are calculated by performing a linear 
interpolation between the appropriate solar flux periods 


reported by Torr et al. (1979) based on the value of F10.7 


(x107 24 Rp ane Hz ~). Thus the EUV flux is scaled according 
Or arene re Since only one measurement of X-ray flux is 
available, interpolation 1s not possible. However, in order 


to explain the observed solar cycle variations of NO, Barth 
et al. (1988) suggested that the soft X-ray flux exhibits a 
larger variation with solar activity than EUV. For “time 
reason, Siskind added a multiplicative X-ray flux scaling 
factor, FRACXR, to the PEGFAC code. The X-ray fluxes used 
by the model are scaled by the multiplicative factor FRAG 
It is the goal of this thesis to determine the relationship 
between FRACKR and F10./7. 

PEGFAC uses the flux data to calculate the 
photoionization rates of N,, 0+, and O, and the _ photo- 
dissociation rates of N, for each model "hour" of every day 


for which the model is run. As with MSIS-83, PEGFAC stores 


the results as 183 element column vectors. For each "hour", 
four of these vectors are generated. The first  thimee 
contain the photoionization rates’ for N5, 05, and oO, 
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respectively. The fourth contains the N. photodissociation 
rates. tise eCdch hour is “asi683) x 4 array. Each day 
produces 24 such arrays. 
c. Direct Inputs to the Model 
The neutral atmosphere, photoionization rates 


for N and O, and the photodissociation rates for No 


eo 
are calculated exterior to the model. The data files 
created by these two procedures are used as inputs to the 
model program. The density profiles for the minor species 
are set to a generic initial condition, and the model is 
allowed to run in a time independent mode for a sufficient 
length of time to allow this initial minor atmosphere to 
reach the steady state conditions which correspond to the 
Starting date of the run. The model "solves" the continuity 
equation for each of the species involved in the nitric 
oxide chemistry every 15 minutes using the numerical methods 
previously described. When the sun is within 15 degrees of 
the horizon the atmosphere is rapidly changing, therefore 
the time increment of the model switches to five minutes. 
As with MSIS-83 and PEGFAC, the model stores the density 
profile for each minor constituent in a 183 element column 
vector. 
3. Outputs From the Model 

The program is written to calculate the density 

profile for the minor species at 15 minute intervals, or at 


five minute intervals if the sun is within 15 degrees of the 


Zo 


hor? zone For the purpose of this thesis, the output has 
been structured to extract the nitric oxide density at 110 
km for each day at 1500 local time. This is because the SME 
satellite data was taken at an altitude of 110 km at 1500 
local time each day. The output is stored in a vector 
containing N elements, where N is the number of days for 


which the model was run. 
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V. ANALYSIS 


Chapters III and IV introduced the theory and modeling 
technique that is used in this thesis. This chapter will 
present the experimental approach used to determine the 


relationship between FRACXR and F10.7. 


A. SME OBSERVATIONS 

Data from the Solar Mesosphere Explorer satellite was 
provided by Dr. Charles A. Barth of the Laboratory for 
Atmospheric and Space Physics, University of Colorado 
(telephone conversation 17 November, 1988). Equatorial 
latitude nitric oxide data was used for this thesis. This 
data is the result of averaging measurements from _ two 
satellite orbits between 20° north latitude and 20° south 
latitude. The first orbit took measurements from 20°-60° 
east longitude, while the second were between 105°-147° west 
longitude. All measurements were taken at an altitude of 
110 km at 1500 local time. The 110 km nitric oxide density 
for the period January 6, 1982 to April 19, 1985 is shown in 
megure 5.1. 

1. Analysis of Original Model 

Figure 5.2 shows a comparison between SME data and 

model calculations prior to the addition of Siskind’s FRACXR 
parameter. It clearly shows that the model was inadequate 


in predicting both the absolute magnitude and variations in 


on 


NO Density (x10" cm”) 
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Figure 5.1. Nitric oxide density at 110 km. 

SME measurements from January 6,1982 to 

ADE1 1 197 965, 
the nitric oxide density which were observed by the SME 
Satellite. The reason for this inadequacy can be explained 
by taking a closer look at the relevant photochemistry. 
Originally, an increase in the EUV flux would increase the 
rate at which the production of N(“D) occurs (see R12 and 
R27 in Appendix). However, this increase in EUV flux would 
also cause one of the principal loss reactions for nitric 
oxide (R22) to increase its rate due to an increase in the 
ionized molecular oxygen density. The increase in n( 7D) 
(and hence NO) was thereby partially offset by the increased 


loss of NO via R22. A process was sought which would allow 
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Figure 5.2. Comparison between SME observations 

and the photochemical model prior to FRACKXR. 
the model to increase the production of N(“D) without 
causing a corresponding decrease in NO density via the 
loss reaction R22. This led to the modification of the 
PEGFAC code by Siskind as discussed in Chapter IV. 

The reason behind the incorporation of the 
multiplicative X-ray scaling factor is now clear. Soft 
X-rays which penetrate to the E-region photoionize the major 
constituents and produce energetic photoelectrons = as 
discussed in Chapter III. The increase in photoelectrons 
drives the production of N(*D) via R12 and R17, but does not 


affect R22. In order for this mechanism to work, the 


ft 


33 


variation in solar X-radiation must be greater than the 


variation in solar EUV radiation. (Barth et al., 1988) 


B. DETERMINATION OF FRACXR 

Initial estimates of the magnitude of FRACXR’ were 
provided by Siskind (telephone conversation November 9, 
1988). He suggested that FRACXR take on values near one for 
conditions of solar minimum and 50 for solar maximum. This 
provided a starting point for the determination of the 
relationship between FRACXR and solar cycle. 

1. Assumptions 

Two major assumptions were made for the purposes of 
simplifying the data analysis and reducing computer time. 
The first assumption is that auroral electrons do not 
influence the nitric oxide chemistry. This assumption is 
based on the fact that all computations and observations 
were at equatorial latitudes. Tt 1s justified by the work 
of Siskind (1988), who showed that for equatorial latitudes 
there is little enhancement of NO density due to an auroral 
storm. 

The second assumption was that the computation of 
nitric oxide density would be independent of the longitude 
at which the model was run. This allowed the model to be 
run at only one longitude for the computation of each NO 
density time sequence, instead of averaging two sets of 
measurements as was done for the SME data. To justify this 


assumption, the model was run at both 40 degrees east 
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longitude and 126 degrees west longitude for the same five 
day period in 1982. Figure 5.3 shows the result of this 
procedure. The two cases deviated from each other by less 
than 0.2%. Based on this result, all subsequent model runs 


were done only at 40 degrees east longitude. 
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Figure 5.3. Model [NO] calculations for 


two different longitudes. 
2. Experimental Approach 
The SME data was searched to locate periods where 
the value of F10.7 remained constant for at least three 
days. Four such periods were initially found corresponding 
to F10.7 (x10~°2 wm? uz~+) values of 70, 90, 128, and 201. 


The model was then run with various values of FRACXR for 
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each of the four periods. Figures 5.4 through 5.7 show the 
results of these runs. This initial data is summarized in 
Table II and shown graphically in Figure 5.8. Uncertainties 
are estimated from the data directly without measuring them 
explicitly since only one set of calculations is possible 


each time the model is run (see e.g., Bevington,1969). 


TABLE II 
SUMMARY OF INITIAL FRACXR DATA 


PeeOre) 


(olen (Satna 





It can be seen from Figure 5.8 that the first three 
data points are nearly linear, while the point corresponding 
to an F10.7 value of 201 lies significantly lower than would 
be expected uSing a linear relationship based on the first 
three points. To resolve this uncertainty, a fifth data 
point was taken for an F10.7 value of 242. =Although this 
value was constant for only two days, the FRACXR value at 
this point would determine if a linear relationship was 


probable. The result of this run is shown in Figure 5.9. 
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Figure 5.4. Model [NO] densities for F10.7 = 7/0. 
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Figure 5.5. Model [NO] densities for F10.7 = 90. 
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Figure 5.7. 
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Figure 5.8. Initial FRACXR Data. 
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Figure 5.9. Model [NO] densities for F10.7 = 242 
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3. Fitting the Data 
The addition of the fifth data point at F10.7 = 242 
is reflected in Table III, a summary of the final FRACXR 
data. The problem now is to find a function (or combination 
of functions) which best fits the data. This function can 
then be used to determine the value of FRACXR on a daily 


basis for input to the model. 


TABLE TIS 
SUMMARY OF FINAL FRACKR DATA 


Fors 7 


(x107 22 Wm FRACXR 





Several fitting functions were tested. The criteria 
for determining the best fit was that function which 
minimized x? and was monotonically increasing in FRACKR. 


2 


"The definition of x as the ratio of the estimated variance 


5? to the parent variance o*...makes it a convenient measure 
of the goodness of fit." (Bevington, 1969). A weighted 
least-Squares fit to a third degree polynomial is shown in 


Figure 5.10, and a weighted least-squares quadratic fit to 
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the first four points is shown in Figure 5.11. Weighting 


was determined for each data point by 





where o; was taken to be the estimated uncertainty. Figure 
5.12 shows a least-squares logarithmic fit, and Figure 5.13 
shows a combination of two linear least-squares fits. 


Results are summarized in Table IV. 


TABLE IV 
SUMMARY OF FITTING FUNCTIONS 


FUNCTION 


THIRD DEGREE POLYNOMIAL 


QUADRATIC 


LOGARITHMIC 


DOUBLE LINEAR FIRST LEG 
SECOND LEG 
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FRACXR = A + Bx(F10.7) + Cx(F10.7)" + Dx(F10.7)3 


A = -—-48.497 

B = 0.98430 

C =>-4.42519E—03 
= D = '7.13504E-08 
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Figure 5.10. Weighted Third Degree Polynomial Fit. 


FRACKR = A + Bx(F10.7) + Cx(F10.7)* 


A = —41.216 
B = 0.71832 
C =>-1.95636E-—-03 
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Figure 5.11. Weighted Quadratic Fit. 
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Figure 5.12. Weighted Logarithmic Fit. 


FRACKXR = A + Bx(F10.7) 
A= 14.930 
B= 4.82255E—02 


FRACKXR = A + Bx(F10.7) 
A= —24.513 
B= 0.35140 
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Figure 5.13. Weighted Double Linear Fit. 
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4. Statement of Empirical Relationship 
Based on the results of the analysis presented in 
the preceding section, the empirical formula for the scaling 
factor FRACXR which relates X-ray flux and EUV flux in the 


E-region of the ionosphere is 


FRACXR = 14.930 + 0.048x(F10.7) 
for F10.7 > 130 and 


FRACKR = -24.513 + 0.351x(F10.7) 


for F10.7 < 130, where F10.7 has units 1077% wm? H2--. 
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VI. RESULTS 


In Chapter V, an empirical expression for FRACXR was 
proposed. In this chapter, the expression is incorporated 
into the one dimensional diffusive photochemical model to 


compute nitric oxide for two 60 day time periods. 


A. COMPARISON OF MODEL AND SME DATA 

Two 60 day periods were chosen to further test the 
model. Each period was chosen to correspond to a time where 
the SME data showed distinct periodic variations in nitric 
oxide density. The first 60 day period was from 82213 to 


-22 W m2 


82273, and was chosen because the F10.7 value (x10 
Hz 1) ranged from a minimum of 120 to a maximum of 214 and 
thus would correspond to FRACXR values taken from the upper 
curve of Figure 5.13. The result of this 60 day run as well 
as the daily value for F10.7 is shown in Figure 6.1. Model 
calculations without uSing X-ray scaling are included for 
comparison. The second 60 day period was from 84210 to 


22 W m2 


84180 , and was chosen because the F10.7 value (x10 
ie) ranged from a minimum of 90 to a maximum of 183 and 
thus would correspond to FRACXR values taken from the lower 
curve of Figure 5.13. The result of this 60 day run, along 


with the F10.7 values and unscaled model calculations, is 


shown in Figure 6.2. 
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Figure 6.1. Comparison of SME Data to Model 
for Period of Solar Maximum. 
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Figure 6.2. Comparison of SME Data to Model 
for Period of Solar Minimun. 


47 


1. Discussion of Results 

The average daily deviation of the model calculat- 
ions for nitric oxide density from the SME satellite 
measurements was 13.92% for the first 60 day run and 24.27% 
for the second 60 day run. This is not surprising, since 
the first run was based on a FRACXKR relationship which has a 
38% probability of being the correct “uneeron (x? « 2759 
while the second run was based primarily on a eFRACKXR 
relationship which only has a 4.8% probability of being the 


correct function (x? = 3.95). Figures 6.1 and 6.2 show that 


variations of X-ray flux by a factor of about 30 over ene 
solar cycle can explain the observed variation in) jmijege 
oxide density. 


Calculated nitric oxide densities for both 60 day 
runs are plotted against F10.7 to test the correlation 
between the two data sets in Figure 6.3. “An inspection of 
the figure shows that correlation does exist and follows the 
same "double linear" type of relation as exists between 
FRACXR and F10.7. A least-Squares fit to each branch of 
this scatter plot yields lines with regression coefficients 
of 0.55 and 0.90 for the upper (F10.7 greater than 130) and 
lower (F10.7 less than 130) parts, respectively. This 
strongly suggests that FRACXR is incorrectly calculated for 
values of F10.7 greater than 130. 

An inspection of Figures 6.1 and 6.2 reveals that 


the derived relationship between FRACXR and F10.7 does not 
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Figure 6.3. Correlation of Calculated NO 
Densities and F10.7. 


the model to correctly calculate NO densities for values of 


-22 2 -1 


Bio.7 (x0 Wm" Hz =) £f£rom about 140 to 160. The 


computed NO densities for these values are too high due to 


the artificially high production rate of N(* 


D) from the 
increased X-ray flux. This suggests that a smoother 
transition exists between the lower and upper curves of 


Figure 5.13 than proposed by the "double linear" fit. 
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VII. CONCLUSION 


A. SUMMARY 

A theoretical and observational study of the production 
of nitric oxide in the E-region of the ionosphere has been 
presented. The observational basis was the NO measurements 
made by the Solar Mesosphere Explorer. satellite. An 
empirical expression for the variation in soft solar x-ray 
flux as a function of the 10.7 cm solar flux was determined. 
This variation was in the form of a multiplicative scaling 


factor, FRACXR, between X-ray and EUV flux given by 


FRACXR = 14.930 + 0.048x(F10.7) 
for F10.7 greater than 130 and 


FRACXR = -24.513 + 0.351x(F10.7) 


for F10.7 less than or equal to 130, where F10.7 has units 
22 Sin See The 10.7 cm flux was chosen as_ the 
independent parameter because it is used in the neutral 
atmosphere model as an indicator of solar EUV radiation. 

The expression which calculates FRACXR was incorporated 
into a one dimensional diffusive photochemical model to 
compute nitric oxide densities. Model calculations for NO 
densities showed significant improvement in both magnitude 


and variation. These model calculations were compared to 


observational measurements from the SME satellite. Results 


=) 9, 


of this comparison showed that model calculations using the 
above X-ray scaling factor were in much better agreement 
with the observed SME data than the calculations made 
without X-ray scaling. The conclusion from this invest- 
igation is that variations of soft solar X-ray flux by a 
factor of 30 over the solar cycle can explain the observed 


variation in nitric oxide densities. 


B. SUGGESTIONS FOR FURTHER RESEARCH 

Figure 7.1 shows a scatter plot of observed nitric oxide 
densities versus 10.7 cm solar flux. Barth et al. (1988) 
pointed out that the observed nitric oxide densities show a 
strong correlation with Lyman-a flux over solar cycle. This 
is shown in Figure 7.2. Examination of these two figures 
suggests a stronger correlation of NO density to Lyman-a 
than to F10.7. It is suggested that further investigations 
be made into model calculations of NO density based on an 


X-ray scaling factor which is a function of Lyman-ca. 
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Figure 7.1. Correlation of NO and F10.7. 
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Figure 7.2. Correlation of NO and Lyman-ca. 
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APPENDIX 


NITRIC OXIDE REACTIONS IN THE LOWER THERMOSPHERE 


This appendix lists all of the reactions used in the 
Cleary one dimensional photochemical model mU965). 
Reactions between neutral species are listed first, followed 
by 1on-neutral reactions. Equation numbers’ shown here 
reflect the scheme used in Cleary (1986), except for R39 


which was added later. 


A. NEUTRAL SPECIES 


1. Photoionization 


O+hv 30 +e Rl 
0, + hv > On. + e R2 
N, + hv » Ne + e R3 
O+ hv » a4 + e R4 
n(4s) + hv 3 NT + RB 
0, +hv?0+0 +e RY 
N5 + hv > n(4s) +N +e R10 
NO + hv + NO’ +e R13 


2. Photodissociation 


N, + hv 4 n(4s) + N(4s) R14 
N, + hv 4 N(“D) + N(2D) R15 
NO + hv 3% O + n(4s) R16 


23 


3. Radiative Relaxation 


nN(“D) 3 N(4S) + hv R32 
4. Quenching 

N(“D) + e tN(‘S) +e R30 

nN(“D) + 0 » N(4S) +0 R31 


5. Bimolecular Reaction 


n(4s) +0 


5 > NO + O R28 
N(“D) + 0, # NO +0 R29 
N(7P) aa 0, > NO + O R39 
N(“S) + NO % Ny +0 R33 
N(“D) + NO >» No + O R34 
6. Energy Transfer or Photoelectron Impact 

* + 
O+e %F%F0 + 2e R5 

* ere 6 
05 +e >» 0-5 + 2e R 

; eee, R7 
N5 +e 9% N, + 2e 

* 4 4 . 
N, +e 3% N( SS) + N( S) + e Ril 

* 2 Zz 
No + e 2% N('D) + N(“D) + e R12 

B. IONS 


1. Charge Transfer 


a +. 

No + 0.5 = 0. + N, R18 
+ “ 

O + 0, => 0. + O R21 
aa + 

0. + NO» NO + 0, R22 

0(“p)* + No => No” + O R35 

O( “Di eo» aC manG R36 
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Ion-Atom Interchange 


a 
Z 


Owns N, > not + n(4s) 
+ 
Z 


nt + QO 


w.t + 0 3 Not + N(“D) 


O + n(4s) > No* + O 


>? Not + 0 


Dissociative Recombination 


+ 
2 


One hee io) aero 


N +eo.> N( 7D) + N(“D) 


not +e4%0+ N(4s) 


No* +er? 0 + N( 7D) 
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R17 
R20 
R23 


R25 


R19 
R24 
R26 


R27 
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